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Abstract 

We discuss the nature of quasi-stationary states (QSS) with non-Boltzmannian dis- 
tribution in systems with long-range interactions in relation with a process of incom- 
plete violent relaxation based on the Vlasov equation. We discuss several attempts 
to characterize these QSS. We show that their distribution is non-universal and 
explain why their prediction is difficult in general. 

Key words: Vlasov equation, long-range interactions, quasi-stationary states 



1 Quasi-stationary states: a generalized thermodynamics? 



It has been observed in many domains of physics [1] that Hamiltonian systems 
with long-range interactions spontaneously organize into coherent structures 
which persist for a long time. Some examples are galaxies in astrophysics, jets 
and vortices in 2D geophysical flows, clusters in the HMF model etc. These 
quasi- stationary states (QSS) are usually not described by the Boltzmann 
distribution. To account for this striking observational fact, some authors have 
proposed to replace the Boltzmann entropy by the Tsallis entropy 

S q [f} = -^- l ju q -f)dvdv. (1) 

The reason advocated is that the system is non-extensive so that standard 
thermodynamics may not be applicable [2]. The maximization of the Tsallis 
entropy at fixed mass and energy leads to g-distributions of the form f(r, v) = 
[fj, — (3(q — l)e/g]9 rrT , where e = v 2 /2 + $(r) is the individual energy and /x, 
(3 are Lagrange multipliers. There are situations where these distributions 
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provide a good fit [2,3] of the QSS. However, there exists other situations [4] 
that are described neither by the Boltzmann nor by the Tsallis distribution 
[5]. We want to show that the prediction of the QSS is relatively complicated 
and explain why. To that purpose, we use classical methods of kinetic theory 
relying on the Vlasov equation [1]. 



2 Kinetic theory: the importance of the Vlasov equation 



To understand the physics of the problem, we first have to develop a ki- 
netic theory of systems with long-range interactions [1,6]. Such kinetic the- 
ories have been developed for stellar systems, two-dimensional vortices, the 
HMF model etc. They usually lead to a kinetic equation for the distribution 
function /(r, v, t) of the form 

df _df df df 1 

^ = ^ + V -^- V$ -^ = A^ W) - (2) 



The left hand side, called the Vlasov term, describes an advection in phase 
space due to the mean-field potential $(r, t) = J u(\r — r'\)f(r',v',t)dv'dr' 
where u is a binary potential of interaction. The right hand side takes into 
account the effect of "collisions" (more generally correlations) between parti- 
cles and depends on the number N of particles. For long-range interactions, 
it usually scales like N~ 5 with S > 1 (in general 5=1 but for ID systems 
5 > 1). This term, due to finite N effects (graininess) , is responsible for the 
collisional relaxation of the system towards ordinary statistical equilibrium. 
Indeed, in general, the collision term cancels out only for the Boltzmann dis- 
tribution: Q(f e ) = <-> f e = Ae~ l3me . However, due to its dependence with 
the number of particles, the collisional relaxation time is of order tR ~ N 5 to 
where to is the dynamical time. In general, this timescale is huge and does not 
represent the regime of most physical interest. In particular, the QSS men- 
tioned previously form on a timescale t ~ to <C tR for which the collision 
term can be neglected in a first approximation. In that case, the evolution 
is collisionless and described by the Vlasov equation: df /alt = 0. Now, it has 
been understood, first in astrophysics by Henon, King and Lynden-Bell in the 
1960's, that the Vlasov equation, when coupled to a long-range force like the 
gravitational force (Vlasov-Poisson system) can undergo a form of collision- 
less relaxation on a very short timescale ~ to- This is called violent relaxation 
[7]. This general process explains the ubiquity of long-lived QSS in systems 
with long-range interactions. The fine-grained DF /(r,v,t) which is solution 
of the Vlasov equation never achieves equilibrium but develops intermingled 
filaments at smaller and smaller scales due to phase mixing. However, if we 
locally average over these filaments, the coarse-grained DF /(r, v,t) is ex- 
pected to converge toward a steady state which is a stable stationary solution 
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of the Vlasov equation. This is called weak convergence in mathematics. The 
kinetic theory explains that the lifetime of the QSS scales as a power N s 
of the number of particles. The kinetic theory also explains that when the 
t — > +00 limit is taken before the iV — > +00 limit we get the Boltzmann sta- 
tistical equilibrium state (collisional relaxation) but when the N — > +00 limit 
is taken before the t — > +00 limit we get a QSS which is a stable stationary 
solution of the Vlasov equation usually different from the Boltzmann distri- 
bution (collisionless regime). This non-commutation of the limits N — > +00 
and t — > +00 has been observed numerically by [2] although they did not 
give an explanation in terms of the Vlasov equation, as was done in [1] (see 
Sec. 6.5) but rather in terms of Tsallis generalized thermodynamics. On short 
timescales ~ to, the system undergoes violent relaxation and reaches a sta- 
tionary solution of the Vlasov equation (on the coarse-grained scale) . Since the 
Vlasov equation admits an infinite number of stationary solutions, the solution 
(QSS) effectively selected by the evolution is difficult to predict (see below). 
On intermediate timescales to < t < £r, the system passes by a succession of 
quasi- stationary states that are quasi-stationary solutions of the Vlasov equa- 
tion f(e,t) slowly evolving with time due to "collisions" (finite N effects). In 
astrophysics, this slow collisional evolution is governed by the orbit-averaged- 
Fokker-Planck equation. On a longer timescale ~ t R , collisions finally select 
the Maxwell distribution among all stationary solutions of the Vlasov equa- 
tion (the fate of gravitational systems is peculiar due to the absence of strict 
statistical equilibrium and the process of evaporation or collapse). These dif- 
ferent regimes have been observed for different physical systems [1,6] and are 
illustrated numerically by Yamaguchi et at. [8] for the HMF model. 



3 Lynden-Bell's theory of violent relaxation 

In a seminal paper, Lynden-Bell [7] tried to predict the QSS resulting from 
violent relaxation assuming that the system mixes well (ergodicity) and us- 
ing arguments of statistical mechanics. However, the statistical mechanics of 
the Vlasov equation is peculiar because of the presence of an infinite set of 
constraints: the Casimirs Ih — I h(f)drdv (for any function h) which contain 
all the moments I n — J f n drdv of the distribution function. These are sort of 
"hidden constraints" [5] because they are not accessible from the (observed) 
coarse-grained DF / since / f n drdv 7^ J f n drdv for n 7^ 1. Therefore, the 
proper density probability to consider in the statistical theory of violent re- 
laxation is p(r, v, rj) which gives the density probability of finding the level of 
DF / = 77 in (r, v) in phase space. Note that we make the statistical mechanics 
of a field, the distribution function, not the statistical mechanics of discrete 
particles. From this density probability, we can construct all the coarse-grained 
moments f n = J prj n drj including the coarse-grained DF f — J prjdr). We can 
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now introduce a mixing entropy from a combinatorial analysis [7,5] like in 
Boltzmann's traditional approach. In the context of the Vlasov equation, this 
entropy is a functional of p(r, v, r/) of the form 



Sl.b.Ip] = ~ p(r, v,rj) hi p(r,v,rj)drdvdr] 



(3) 



The Lynden-Bell entropy (3) is the proper form of Boltzmann entropy taking 
into account the specificities of the Vlasov equation. Assuming ergodicity (effi- 
cient mixing), the QSS is obtained by maximizing Sl.b. [p] at fixed mass, energy 



where a and (3 are the usual Lagrange multipliers associated with the con- 
servation of mass M and energy E while = ex P{ _ En a n J ] n } takes into 
account the conservation of all the Casimirs invariants I n [5]. We note that, 
even if the system mixes well (ergodicity), the coarse-grained DF predicted 
by Lynden-Bell may differ from the Boltzmann distribution: /l.b. ( e ) Ae^^ 
(in general) due to the presence of the Casimir invariants. In the simplest case 
(two-levels approximation) the DF predicted by Lynden-Bell is similar to the 
Fermi-Dirac statistics: ff.B. = ?7o/(l + e^ - ^) [9]. More generally, the coarse- 
grained DF (4) is a sort of superstatistics [5,10] as it is expressed as a superposi- 
tion of Boltzmann distributions (universal) weighted by a non-universal factor 
x(v) depending on the initial conditions. Furthermore, like for the Beck-Cohen 
superstatistics, the coarse-grained DF (4) maximizes a "generalized entropy" 
in /-space S[J] = -JC(f)drdv with C(f) = - f J [{\xix)'}-\-x)dx (where 
x($) = /^xfl)^''*^) at fixed mass and energy [5]. The distribution (4) 
and the corresponding entropy S[f] are non-universal. However, a general pre- 
diction of the Lynden-Bell statistical theory [7] is that the QSS is a stationary 
solution of the Vlasov equation of the form f = f(e) with / (e) < 0: the DF 
depends only on the energy and is monotonically decreasing. Furthermore, 
/( r ) v ) < max riV /(r, v, t = 0): the coarse-grained DF is bounded by the max- 
imum value of the initial DF. 



4 Incomplete violent relaxation 

The Lynden-Bell DF (4) is the proper prediction of the QSS when mixing is 
efficient (ergodic) during violent relaxation. There are situations where the 
Lynden-Bell prediction works relatively well, e.g. [11]. However, it has been 
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.^( 7 y) e - r ?(/ 3e + a ) 
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recognized in many other occasions [4] that mixing is not efficient enough to 
sustain the hypothesis of ergodicity on which the theory is built so that the 
prediction of Lynden-Bell fails in practice: fcjss ( m general) [1]. 

This is particularly obvious in the case of self-gravitating systems because the 
Lynden-Bell entropy has no maximum at fixed (finite) mass and energy [7,9]. 
What can we do to account for incomplete relaxation! 

A first possibility would be to change the form of entropy. For example, we 
could try to use the "generalized thermodynamics" of Tsallis. However, as 
explained previously, the proper density probability is p(r, v, r/) so that the 
proper form of Tsallis entropy in the context of Vlasov systems is [5] : 

S q [p] = l — I [p q (r, v, 77) - p(r, v, 77)] drdvdi], (5) 



instead of (1). For q — 1 it returns the Lynden-Bell entropy (3). In this line 
of thought, q would measure the efficiency of mixing (q — 1 if the evolution is 
ergodic). For q ^ 1, the functional (5) could describe non-ergodic behaviours 
(incomplete mixing). However, it is not clear why all non-ergodic behaviours 
could be described by a simple functional such as (5). Tsallis entropy may 
describe a certain type of non-ergodicity (fractal or multifractal) with phase- 
space structures but probably not all of them. In particular, observations of 
elliptical galaxies in astrophysics [4] do not favour Tsallis form of entropy since 
elliptical galaxies are not stellar polytropes that would be the prediction based 
on the g-entropy [5]. Therefore, /qss 7^ f q (c) (in general). 

Another possibility is to keep the Lynden-Bell entropy as the most funda- 
mental entropy of the problem but develop a dynamical theory of violent 
relaxation. Indeed, if relaxation is incomplete, we must understand why. Qual- 
itatively, the collisionless relaxation is driven by the fluctuations of the field 
$(r, t). Now, these fluctuations can vanish before the system had time to relax 
completely so that the system can remain frozen in a stationary state of the 
Vlasov equation which is not the most mixed state. By using a phenomeno- 
logical Maximum Entropy Production Principle (MEPP), we have proposed 
in [12] to describe the out-of-equilibrium evolution of the probability density 
p(r, v, i], t) by a relaxation equation of the form 



If the diffusion coefficient were constant, this equation would relax towards 
the Lynden-Bell distribution. However, it is argued in [12] that D is not con- 
stant. Indeed, the relaxation is driven by the fluctuations of the field $, itself 
induced by the fluctuations of /, so that the diffusion coefficient should be pro- 
portional to p = f 2 — f and vanish in the regions of phase space where these 



5 



fluctuations vanish. Moreover, as the system approaches (quasi)-equilibrium, 
the fluctuations of the field are less and less efficient so that the diffusion 
coefficient should also decay with time. For these reasons, it can become very 
small -D(r, v, t) — > in certain regions of phase space (where mixing is not very 
efficient) and for large times (as the fluctuations weaken). The diffusion coef- 
ficient can also rapidly decay with the velocity. The vanishing of the diffusion 
coefficient can "freeze" the system in a subdomain of phase space (bubble) and 
account for incomplete relaxation and non-ergodicity. The relaxation equation 
(6) should then tend to a distribution which is only partially mixed and which 
is usually different from the Lynden-Bell and the Tsallis distributions. How- 
ever, this approach demands to solve a dynamical equation (6) -smoother than 
the Vlasov equation- to predict the metaequilibrium state. The idea is that, 
in case of incomplete relaxation (non-ergodicity), the prediction of the QSS 
is impossible without considering the dynamics: it depends on the "route to 
equilibrium" . 

Maybe, we have to accept that, in general, the QSS is unpredictable in case 
of incomplete relaxation. We can expect, however, that /qss(i", v) i s a s ^ a_ 
ble stationary solution of the Vlasov equation. We are thus led to construct 
stable stationary solutions of the Vlasov equation is order to reproduce obser- 
vations. The Vlasov equation admits an infinite number of stationary solutions 
given by the Jeans theorem, but not all of them are stable. Of course, only 
stable solutions must be considered and their selection is a difficult prob- 
lem. If we restrict ourselves to DF of the form / = /(e) with /'(e) < 
depending only on the energy (in astrophysics, such distributions character- 
ize a subclass of spherical stellar systems), it is possible to provide a simple 
criterion of nonlinear dynamical stability. Such DF extremize a functional of 
the form H[f] = — J C(f)drdv, where C is convex, at fixed mass and en- 
ergy [13]. Indeed, the first variations SH — (36 E — aSM = yield a DF of 
the form / = F((3e + a) with F(x) = {C'Y l (—x) montonically decreasing. 
If, furthermore, the DF maximizes this functional (at fixed E, M), then it is 
nonlinearly dynamically stable with respect to the Vlasov equation [13,8,14]. 
The intrinsic reason is that H[f], E[f] and M[f] are individually conserved 
by the Vlasov equation. Therefore, if / (r, v) is the maximum of H[f] (at 
fixed E, M), a small perturbation /(r, v,t) will remain close (in some norm) 
to this maximum. It is important to note that this criterion of nonlinear dy- 
namical stability is remarkably consistent with the phenomenology [5] of vio- 
lent relaxation if we view the relevant DF as the coarse-grained DF. Indeed, 
during mixing df/dt ^ and the functionals H[f] = — J C(f)drdv calcu- 
lated with the coarse-grained DF increase (-H decrease) in the sense that 
H[f(r, v, t)] > H[f(r, v, 0)] for t > where it is assumed that initially the 
system is not mixed: /(r, v, 0) = /(r, v, 0). Because of this property similar 
to the Boltzmann if -theorem in kinetic theory 1 , H[f] are called generalized 



Note that the Vlasov equation does not single out a unique H-function contrary 
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iJ-functions [13]. By constrast, E[f] and M[f] are approximately conserved. 
Therefore, this generalized selective decay principle [5] (decrease of —H[f] 
at fixed E, M) due to phase mixing and coarse-graining can explain how f 
can possibly reach a maximum of H at fixed mass and energy (while H[f] is 
rigorously conserved on the fine-grained scale). After mixing df/dt = and 
the functionals H[f] (as well as E[f] and M[/]) are conserved by the coarse- 
grained flow. Therefore, if / has reached (as a result of mixing) a maximum 
of H, it will be nonlinearly dynamically stable with respect to coarse-grained 
perturbations (after mixing) in virtue of the above-mentioned stability result. 

In reality, the problem is more complicated because the system can converge 
toward a stationary solution of the Vlasov equation which does not depend on 
the energy e alone, and thus which does not maximize an if-function at fixed 
E, M. For example, the velocity distribution of stars in elliptical galaxies is 
anisotropic and depends on the angular momentum J = |r x v| in addition 
to energy e. Furthermore, real stellar systems are in general not spherically 
symmetric so their DF does not only depend on e and J. Therefore, more 
general stationary solutions of the Vlasov equation must be constructed in 
consistency with the Jeans theorem. Elliptical galaxies are well relaxed (in the 
sense of Lynden-Bell) in their inner region (leading to an isotropic isothermal 
core 2 with density profile ~ r~ 2 ) while they possess radially anisotropic 
envelopes (with density profile ~ r -4 ). Stiavelli & Bertin [4] introduced an /(°°) 
model of the form /(°°) = A(—e) 3 / 2 e~ ae ~ cJ I 2 for e < (and / = otherwise) 
based on the possibility that the a priori probabilities of microstates are not 
equal due to kinetic constraints. This model reproduces many properties of 
ellipticals but it has the undesired feature of being "too isotropic". Then, 
they introduced another model based on a modification of the Lynden-Bell 
statistical theory. They considered the maximization of the Boltzmann entropy 
(in Lynden-Bell's sense) at fixed mass, energy and a third global quantity 
Q = J J u \€\~ 3u / 4 fdrdv which is argued to be approximately conserved during 
violent relaxation. This variational principle results in a family of models 

f( u ) = A exp[— at — d (</ 2 /|e| 3 / 2 ) ]. These models are able to fit products of 
TV-body simulations over nine orders of magnitude in density and to reproduce 
the de Vaucouleur's R 1 ^ law (or more general R x l n laws) of ellipticals. The 
introduction of additional constraints in the variational principle could be a 
way to take into account effects of incomplete violent relaxation. 

We note finally that the Vlasov equation can have a very complicated, non- 



to the Boltzmann equation (the above inequality is true for all //-functions) and the 
time evolution of the //-functions is not necessarily monotonic (nothing is implied 
concerning the relative values of H(t) and H(t') for t,t' > 0). 

2 One success of Lynden-Bell's theory of violent relaxation is precisely to explain 
the isothermal cores of elliptical galaxies without recourse to "collisions" whose 
effect manifests itself on a much longer timescale. 
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ergodic, dynamics. For example, in the gravitational ID Vlasov-Poisson sys- 
tem, phase-space holes which block the relaxation towards the Lynden-Bell 
distribution (incomplete relaxation) have been observed [15]. In that case, the 
system does not even relax towards a stationary state of the Vlasov equa- 
tion but develops everlasting oscillations. Rapisarda & Pluchino [16] have also 
observed transient phase-space structures in their iV-body simulations of the 
HMF model and they have proposed an interesting analogy with glassy dy- 
namics. These results are not necessarily in contradiction with the Vlasov 
equation (as they beleive), even for the inhomogeneous situations that they 
consider. It would be of interest to check whether these phase-space structures 
can also be obtained by solving the Vlasov equation for the HMF model. 



5 Summary and conclusion 

Non-Boltzmannian distributions appear in the study of Hamiltonian systems 
with long-range interactions. For these systems, the collisional relaxation time 
towards statistical equilibrium (Boltzmann distribution) is huge because it in- 
creases as a power of the number of particles N. Therefore, the evolution of 
the system is described by the Vlasov equation on a very long timescale [1]. 
Now, because of phase mixing and violent relaxation, the Vlasov equation can 
spontaneously lead to the formation of coherent structures: galaxies in astro- 
physics, jets and vortices in hydrodynamics, clusters in the HMF model, bars 
in disk galaxies... These QSS are (nonlinearly) dynamically stable stationary 
solutions of the Vlasov equation which are not necessarily described by the 
Boltzmann distribution. Indeed, the Vlasov equation admits an infinite num- 
ber of stationary solutions and the system can be trapped in one of them. In 
general, the QSS is non-universal as it depends: (i) on the detailed structure of 
the initial conditions (through the Casimirs, in addition to mass and energy) 
and (ii) on the efficiency of mixing (ergodicity). In this Vlasov context, the 
Tsallis distributions are particular stationary solutions of the Vlasov equation 
which correspond to what are called stellar polytropes in astrophysics [5,14]. 
In addition, the Tsallis functional S q [f] is a particular if -function [13], not an 
entropy (which would be a functional S q [p\). Its maximization at fixed mass 
and energy yields a criterion of nonlinear dynamical stability with respect to 
the Vlasov equation [8,13,14], not a criterion of generalized thermodynamical 
stability. A formal thermodynamical analogy [14] can however be developed to 
investigate the nonlinear dynamical stability problem. 
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